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ABSTRACT 

Binary and multiple systems constitute more than half of the tot al stellar popula- 
tion in the Solar neighborhood (|Kiseleva-Eggleton fc Eggletonll200lD. T heir frequent 
occurrence as well as the fact that more than 70 (|Schneider et al.ll201ll ) planets have 
already bee n discovered in such con figurations - most noteably the telluric companion 
of a Cen B (|Dumusque et al.1l2012f) - make them interesting targe ts in the search for 
habitable worlds. Recent studies (jEggl et alll201 2b; Forgan 2012) have shown, that 
despite the variations in gravitational and radiative environment, there are indeed cir- 
cumstellar regions where planets can stay within habitable insolation limits on secular 
dynamical timescales. In this article we provide habitable zones for 19 near S-Type bi- 
nary systems from the Hipparchos and WDS catalogues with semimajor axes between 
1 and 100 AU. Hereby, we accounted for the combined dynamical and radiative influ- 
ence of the second star on the Earth-like planet. Out of the 19 systems presented, 17 
offer dynamically stable habitable zones around at least one component. The 17 poten- 
tially habitable systems contain 5 F, 3 G, 7 K and 16 M class stars. As their proximity 
to the Solar System (d < 31 pc) makes the selected binary stars exquisite targets for 
observational campaigns, we offer estimates on radial velocity, astrometric and transit 
signatures produced by habitable Earth-like planets in eccentric circumstellar orbits. 
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1 INTRODUCTION 

The discovery and confirmation of te r restri a l bodies orbit- 
ins ot h er stars (e.g. Dumuscme et all J2012I) ; iBorucki et alj 
<|2012l ); iBoruckil l|201ll) ; iLeger et all <|2009p ) has generated 
enormous public as well as scientific interest. It has shown 
that after a mere two decades of exoplanetary research find- 
ing potentially habitable worlds around other stars seems 
to be almost within our grasp. Close-by stars and stellar 
systems are thereby premium targets, as they tend to of- 
fer reasonable signal to noise ratios (SNRs) for photome- 
try and radial veloc ity as well as comparatively large astro- 
metric amplitudes jBeauge et al.ll2007l ; iMalbet et al.ll201ll ; 
iGuedes et aT1l200Sl ; lEggl et al.ll2012al ). As more than half of 
the stars in the Solar neighborhood are members of binar y 
or multiple systems (|Kiseleva-Eggleton Sz Eggletonl [200ll ). 
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it is not surprising that more tha n 70 planets in or ar ound 
binary-stars have been discovered ( Schneider et al 1l201ll ) de- 
spite the current observational focus on single star systems. 
Even-though NASA's Kepler mission ha s been quite success - 
ful in finding circumbinary planet s fe.g.lWelsh et all (|2012l ); 
lOrosz et all (|2012l ); iDovle et al] |201ll )) we will focus on 
binary-star systems with potential circumstellar habitable 
zones (HZs) in this study. In fact, most of the planets dis- 
covered in double stars are in these so-called S-Ty pe con- 
figurations l|Rabl fc Dvoraklll988l ; lRoell et al.ll2012l ). where 
the planet orbits one star o nly. The telluric compa nion of a 
Cen B is such an example |Dumusaue et ai]|2012l ). 

An interesting question in this regard is doubtlessly: 
Can S-Type b inary stars harbor habit able worlds? Already 
Huand ((l960f) and iHarringtonl j 19771 ) and more recently 



Forganl ( 2012 ) investigated the effects such configurations 



have on the insola tion hypothetical planets would receive. 
lEggl et all (|2012bl ) (in the following referred to as EG12) 
were able to derive analytic expressions to find HZs in bi- 
nary star systems unifying dynamical and radiative bal- 
ance models for S-Type binary star - planet systems. While 
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the exact manner in which planets form in tight binary 
star systems is still hotly debated in astrophysical litera- 
ture - s ee for instance | Miiller fc Kiev! lf2012h: iBatvgin et alJ 
l|201ll ); [Pa"ardekooper fc Leinhardtl < |201(J >: iThebaultl (|20ilf ) 
and references therein, the discovery of a Cen B b has 
made the existence of terrestrial planets in S-Type binary 
star systems an observational fact. Opinions still differ on 
whether it is theoretically possible that planets in a Cen- 
tauri's HZs can form on stable orbits. Even though clas- 
sical N-body simulations with best case accretion scenar- 
ios seem to be able to produce terrestrial planets near the 
HZs of the a Cent a uri system dQuintana fc Lissauerl l20ld ; 
iGuedes et afll2008h , iThebault et all l|2009l . l2008h concluded 
that even when gas drag is included the encounter veloc- 
ities between kilometer sized planetesimals would lead to 
erosive collisions, thus making constant accretion unlikely. 
However, in their model they did not include a self-consistent 
evolution of the gas disc, nor did they consider planetes- 
imal self gravitation or re-accre tion of collisional debris. 
iPaardekooper fc Leinhardtl (|2010l ) used a self consistent disc 
model with planetesimals. They were able produced ac- 
cretion friendly scenarios when the collision frequency was 
sufficiently high to prevent orbital dephasing. Other pos- 
sible solutions to the problem of high encounter velocities 
range from i ncludi ng planetesimal and embryo migration 
iPavne et all (2009) over mild inclinati on of planetesim al 
discs with respect to the binary's orbit l|Xie et al.ll2010h to 
more realistic radiative m odeling of the system's gaseous 
disc dMiiller fc Klcy 20l|). 



lEggl et al 
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problems in modeling tidal locking of planets in binary sys- 
tems are mentioned in section [7] A summary (section [8} 
concludes this study. 



2 SELECTION OF BINARY STAR SYSTEMS 

We preselected all detached binaries with semimajor axes 
1 < a,b < 100 AU using the stellar orbital parameters pro- 
vided in the WDC, in order to find suitable S-Type sys- 
tems in the Solar neighborhood where Earth-like planets in 
HZs could be detectable. Hereby, we only considered sys- 
tems within a distance of d < 31 pc from o ur Solar Sys- 
tem a s determined by the Hipparchos mission l|van Leeuwenl 
l2007h . Together with the prerequisite that the binaries' or- 
bital elements had to be available, the afore mentioned re- 
strictions reduced the number of admissible double star sys- 
tems to 313. Furthermore, only double star configurations 
with known spectal types of both components were used. 
Peculiar spectra that might have been classified incorrectly 
by Hipparchos like HIP 17544 & 73695 were also excluded, 
narrowing the set of candidate systems from 313 to 35. The 
ultimate selection criterion consisted of calculating the bi- 
naries' periods using bolometric luminosity derived masses 
together with the semimajor axes given in WDS and com- 
paring them to the observed binary periods. The stellar 
bolometric luminosities and masses required for this pur- 
pose were derived as follows: With the distances available 
through Hipparchos data the systems' absolute visual mag- 
nitudes could be calculated. In order to assess the bolometric 
luminosities of the binary sample, we performed bolometric 
corrections (BC) of t he abso l ute v isual magnitudes using 
the polynomial fits bv lFlowerl (|l996h . The required effective 
temperatures were estimated via spectral type and luminos- 
ity using the ATLAS9 ca talog of stellar model atmospheres 
ijCastelli fc Kurucdlioolh . We then calculated the binaries' 
periods using the masses derived via the mass- luminosity 
relations given in ISalaris fc Cassisil (|2005l ). Only those sys- 
tems who's derived periods did not deviate more than 11% 
from the observed periods were selected for the final sample. 
Stellar and orbital parameters for the final set of 19 S-Type 
binary systems are presented in Tabled] 

In the next section we will briefly discuss the main 
points on how to determine HZs for these S-Type binary 
systems. 



show that even if additional Earth- 
like planets in a Centauri do exist, it is not an easy task 
to find them given current observational limitations in ra- 
dial velocity resolution. The RV signal semi amplitude of 
a Cen B b was near the current edge of feasibility with 
S RV ~ 50cm/s , whereas accuracies lower than 10 cm/s 
would be necessary to discover telluric planets in a Cen 
B's HZs. Astrometry is not much more helpful in this case, 
as the necessary astrometric amplitudes to detect habitable 
worlds in a Centauri will on ly be available near th e end of 
the GAIA mission's lifetime l|Hestroffer et al.ll2O10h . 

In this study we tackle the question whether there are 
S-Type systems in the Solar neighborhood that might make 
for easier targets. For this purpose we select 19 S-Type bi- 
nary systems fro m "The Washingto n Visual Double Star 
Catalog" (WDC) (|Mason et alj|2012j ) with well determined 
stellar parameters that lie within a distance of 31 pc from 
the Solar System, and calculate HZs for each s tellar com- 
ponent using the analytic method presented in lEggl et all 
l|2012bh . We provide estimates on the RV and astrometric 
(AM) root mean square (rms) signal strengths expected for 
an Earth-like planet orbiting at the borders of a system's 
HZs. Furthermore, we present likely transit depths for po- 
tentially transiting habitable planets in co-planar S-Type 
double star systems. 

This article is structured as follows: First, we will dis- 
cuss the selection criteria for the 19 systems investigated 
(section [2]). After a brief summary of the main factors that 
determine habitability for terrestrial planets in binary star 
systems (section |3J, the issue of dynamical stability of plan- 
ets in such configurations is addressed in section [4] Our re- 
sults - tables with HZ borders and signal strength estimates 
- are presented and discussed in sections [5] and [6] Current body systems, such as the planet - binary star configurations 



3 HABITABILITY OF EARTH-LIKE PLANETS 
IN S-TYPE BINARY STAR SYSTEMS 

The most pronounced difference between determining clas- 
sical HZs and HZs for Earth-like planets in binary star sys- 
tems lies in the assumption that planetary orbits are ba- 
sically circular . In fa ct, the well known borders defined by 
iKasting et all (|l993l ) are built on the premises that plan- 
etary insolation will change only on stellar evolutionary 
timescales. Thus, the planet is thought to remain more or 
less at the same distance from its host star on a circular or- 
bit. This a ssumption is implicitly made in almost all recent 



bit. 1ms assumption is implicitly made m almos t all recent 
works , e.g. Kane & Gclino (2012); Pier rehumbert fc Gaidosl 



201 if ): Kaltcncggcr & Sassclov (|2011 



). However, in three 
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Table 1. Orbital and stella r parameters of the 19 investigated bi nary star systems. The values of parameters printed in 
bold letters are taken from l|Mason et alj|2012t lvan~L ccuwcn 20q3)i t ne others were derived as described in section [2] The 
binary's eccentricity and semimajor axis are denoted by at, and e\>\ I is the system's inclination to the plane of the sky. A 
binary components' masses are symbolized by Ma and Mg, their respective luminosities by La and Lg and their effective 
temperatures are denoted T £ jja an d T e jjg. Stellar classifications are given in the columns headed "class A" and "class 
B". 



HIP ID 


a b 


eb 


J 


d 


M A 


M B 




La 


L B 


T effA 


T effB 


class A 


class B 


14669 


9.0 


0.14 


96.8 


15.8 


0.56 


0.39 





096 


0.026 


3580 


3370 


M2 


M4 


30920 


4.3 


0.37 


51.8 


4.1 


0.22 


0.08 





007 


0.001 


3370 


3145 


M4V 


M5.5V 


31711 


42.7 


0.34 


93.9 


21.3 


1.03 


0.57 


1 


137 


0.109 


5860 


4060 


G2V 


K7Ve 


44248 


10.4 


0.15 


131.4 


16.1 


1.44 


0.89 


4 


285 


0.638 


6740 


5250 


F3V 


K0V 


45343 


97.2 


0.28 


21.0 


5.8 


0.52 


0.51 





073 


0.067 


3850 


3850 


MOV 


MOV 


51986 


9.9 


0.75 


129.1 


26.8 


1.88 


1.29 


12 


535 


2.790 


6710 


6740 


F4IV 


F3 


58001 


11.7 


0.30 


51.0 


25.5 


2.94 


0.79 


65 


255 


0.397 


9520 


4780 


AOVe 


K2V 


64241 


11.8 


0.50 


90.1 


17.8 


1.30 


1.12 


2 


887 


1.553 


6440 


6360 


F5V 


F6V 


64797 


89.2 


0.12 


93.4 


11.1 


0.73 


0.52 





277 


0.072 


5015 


3715 


K1V 


M1V 


66492 


46.9 


0.61 


36.3 


22.0 


0.59 


0.48 





121 


0.054 


3782 


3647 


M0.5 


Ml. 5 


67422 


32.7 


0.45 


47.4 


13.4 


0.72 


0.65 





273 


0.174 


4560 


4205 


K4V 


K6V 


84425 


7.7 


0.49 


115.2 


30.6 


1.23 


0.86 


2 


267 


0.556 


6280 


5860 


F7V 


G2V 


84720 


91.6 


0.78 


35.6 


8.8 


0.79 


0.50 





393 


0.062 


5570 


3850 


G8V 


MOV 


87895 


2.4 


0.41 


68.0 


28.2 


1.19 


0.90 


2 


031 


0.648 


5860 


4780 


G2V 


K2V 


93825 


32.7 


0.32 


149.6 


17.3 


1.27 


1.25 


2 


570 


2.432 


6200 


6200 


F8V 


F8V 


101916 


15.7 


0.80 


107.0 


30.1 


1.61 


0.37 


6 


794 


0.023 


5745 


4420 


guv 


K2IV 


106972 


5.3 


0.29 


69.4 


24.5 


0.57 


0.43 





105 


0.033 


3370 


3370 


M2 


M4 


114922 


6.7 


0.44 


117.1 


30.8 


0.49 


0.52 





059 


0.073 


3715 


3580 


Ml 


M2 


116132 


42.5 


0.20 


123.5 


6.2 


0.38 


0.20 





025 


0.006 


3370 


3305 


M4 


M5 




[AU] 




[deg] 


[pc] 


[M@] 


[Mq] 


[Lq] 


[L B ] 


[K] 


[K] 







we are investigating, gravitational interactions will alter the 
planetary orbit. 

Perturbation theory of hierarchical triples predicts, that 
the orbit of the inner pair - in our case host-star and 
planet - will experience significant alterations in its eccen- 
tricity, whereas its semimajor axis remain s almost constant 
l|Marchallll99"til ; lGeorgakarakosll2002l . I2003T I . For nearly equi- 
planar systems, the influence of planetary inclination and 
ascending node to the overall dynamics can be considered 
small, they will be neglected in what follows. Even-though 
there may be short periodic variations, some important 
changes in a planet's orbit happen also on secular timescales. 
Secular periods are usually much larger than the planet's 
orbital period. However, they are a lot smaller than stellar 
evolutionary timescales for detached binary systems with 
semimajor axes a;, < 100 AU. It is thus necessary to include 
the effects of changing planetary orbits in our estimates re- 
garding HZs within binary star environments. In their work, 
EG12 confirmed that variations in the planet's orbit are even 
more important for changes in its insolation than the addi- 
tional radiation form the second star! The only exceptions 
to this rule are systems where the second star is much more 
luminous than the planet's host-star (Lb /La > 4, where 
binary component A is the planet's host-star in this case). 
Therefore, a planet's eccentricity is a dominating factor in 
determining habitability. Yet, how eccentric can a planetary 
orbit become, in order to still allow for habitability? 

IWilliams fc Pollard! (|2002T ) concluded that an Earth- 
like atmosphere together with surface oceans can buffer 
the harsh changes between high insolation at periastron 
and long cold phases near apoastron up to eccentricities 
of e p « 0.7, as long as the average insolation is compa- 
rable with the current insolation of the Earth. Although 
planetary eccentricities of such magnitude are usually not 
reached in close S-Type setups (EG12), the region where 



the planet remains within classical insolation boundaries is 
still strongly impacted. In order to distinguish orbital zones 
that are only habitable "on average" and zones where the 
planet will never exceed classical insolation limits, EG12 in- 
troduced three types of HZs for binary star systems: 

Permanently Habitable Zone (PHZ): The PHZ is 

the region where a planet stays within habitable insola- 
tion limits for all times, despite the changes its orbit ex- 
periences due to gravitational interactions with the sec- 
ondary. For this study, we have chosen the classical run- 
away /maximum gre enho use insolation limits as d efined by 
iKasting et ail (|l993h and lUnderwood et al.l (|2003T ) (KHZ) 

Extended Habitable Zone (EHZ): The binary-planet 
configuration is still considered to be habitable when most 
of its orbit remains within the HZ boundaries. This is true if 
the average received insolation plus one standard deviation 
does not put the planet beyond KHZ insolation limits. 

Averaged Habitable Zone (AHZ): Even an elevated 
planetary eccentricity (e < 0.7) may not be prohibitive 
for habitability since th e atmosphere acts as a buffer 
I Williams fc Pollard! [iobl l. if the time averaged insolation 
stays within habitable limits. The AHZ represents such re- 
gions. 

For details on the definition and calculation of PHZ, EHZ & 
AHZ we refer the reader to EG12. We use the interpolation 
formulae given in lUnderwood et al.l (|2003l 1 to calculate ef- 
fective insolation values for the selected stellar types. After 
a brief discussion concerning aspects of dynamical stability, 
the application of the proposed classification scheme to the 
19 selected binary star systems will be presented in the next 
section. 
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4 DYNAMICAL STABILITY OF 

CIRCUMSTELLAR PLANETS IN BINARY 
STARS 

As was briefly mentioned during the introduction, there 
are many open questions regarding the formation of plan- 
ets in double star environments l|Thebaultl l201ll ). How- 
ever, once formed a planet can survive in the dynam- 
ically stable region around one of the binary compo- 
nents - a fact prove n by observed planets in S-Type 
binar y configurations ( Dumusaue et alJ 120121 : iRoell et al.l 
120121 ; iGiuppone etHI l2012h . If the necessary dynami- 
cal prerequisites are fulfilled, even both stars can har- 
bor planets. Generaliz e d dyn a mical investigations such 



_ Holman fe WiegertJ J 19991), IPilat-Lohinger fc Dvorakl 
(2002), semi- analytical JPichardo et al.ll2005l) or analytica l 
approaches [jSzebehelv fc McKenzid 1 19771 : lEggletonl Il983h 
can be used to determine regions where a test-planet can 
remain on a stable orbit on secular dynamical timescales. 
As the setup used in this work consists of a planar bi- 
nary - Earth configuration, the restricted three body ap- 
proach used in the articles mentioned above can be con- 
sidered a reaso nable approximation . We w ill apply the nu- 
merical fit by [ Holman fc WiegertJ (| 19991 ) and results by 
IPilat-Lohinger fc Dvorakl (|2002l ) to find critical semimajor 
axis for circumstcllar motion. 



5 RESULTS 

The different HZs discussed in section [3] are presented for a 
fictitious Earth-like planet in each of the selected double star 
systems (Fig. HI- The region of instability (striped) is also 
marked. The left graph of Fig. [T] represents HZs around the 
primary (S-Type A), and the right graph shows HZs around 
the secondary (S-Type B) (|Whitmire et al.ll998l ). Black (red 
online) denotes regions which are non-habitable due to ex- 
cessive or insufficient insolation, dark gray (yellow online), 
medium gray (green online) and light gray (blue online) rep- 
resent the AHZ, EHZ and PHZ respectively. Dashed and full 
T symbols give the inner and outer bord ers of the classical 
HZ as defined bvlKasting et al.l (Il993l ) and lUnderwood et al.l 
(|2003h (KHZ). lEggl et all |2012bl ) found a good correspon- 
dence between the KHZ and the AHZ, which is also mir- 
rored in the results at hand. Exceptions are the systems 
HIP 58001 & 101916 where the more luminous companion 
shifts the HZs of the less luminous one considerably. Out 
of the 19 selected systems, 17 permit Earth-like planets in 
HZs on dynamically stable orbits around at least one stel- 
lar component. In total, the 17 habitable systems feature 
16 M, 7 K, 3 G and 5 F class stars. Even if the all F and 
M class stars were to be excluded from the list of hosts for 
HZ - either because of their comparatively short lifespans 
ijKasting et al.| [l993) or tidal and radiative effects (see sec- 
tion [7]l - more than 26% of the stars in this sample would be 
capable of sustaining habitable planets on secular dynam- 
ical timescales. If the stars' mass loss via s tellar winds is 
negli gible, and no cataclysmic events occur |Veras fc Tout! 
120121 ). habitability might be given even for stellar evolution- 
ary timescales. 

A detailed listing of HZ-borders as well as expected ra- 
dial velocity (RV) and astrometric (AM) signal strengths 



produced by a terrestrial planet in the selected systems 
is presented in Tables [2] - SJ Maximum and root mean 
squ are (rmsfl signal s trengths have been calculated follow- 
ing [EggTet^aT] (|2012al ). The corresponding equations are re- 
peated in appendix[X]for the reader's convenience. Compar- 
ing AM and RV signal strengths one can see that - current 
observational equipment assumed - RV seems to stand a 
better chance to find Earth-like planets in HZs of nearby 
double stars. With the discovery of a Cen B b the currently 
feasible RV resolution is approximately 50 cm/s. For the de- 
tection of habitable planets in the a Centauri system, how- 
ever, semi-amplitu des around 10 cm/s would be required 
(Eggl ~et al.ll2012al ). Possible candidate systems such as HIP 
14699, 30920, 106972, 114922 or 116132 would offer better 
conditions for finding habitable Earth analogues via RV than 
a Centauri does. 

As 9 out of the 17 potentially habitable systems fea- 
ture M-stars, it is worth mentioning that determining the 
effective insolation a terrestrial planet receives might not be 
enough to claim habitability. In fact, lLammer et all (|201lD 
are convinced that the potentially elevated level of X-ray and 
extreme UV radiation in M-stars might lead to a different 
atmospheric evolution of an Earth-like planet in an M-star's 
HZ, thus preventing the existence of life as we know it. Ulti- 
mately direct observation of the interaction between stellar 
and planetary atmospheres will be necessary to determine 
to which degree planets can remain habitable in the vicin- 
ity of M-typ e stars. The propos ed transit spectroscopy mis- 
sion ECHO l|Tinetti et al.ll2012l ) can be a step in this direc- 
tion, although currently only super-Earths down to 1.5 r® 
around K-F stars are planned to be observed. With RV sig- 
nal amplitudes of ~ 5 — 12 cm/s for potentially habitable 
planets in systems containing Sun-like G stars (HIP 31711 & 
84425), our e stimates are compa rable t o those for a Centaur i 
presented in lEggl et al.1 (|2012al ) and iGuedes et all (|2008h . 
Detecting planets around Sun-like stars would therefore re- 
quire a considerable amount of dedicated observation time 
l|Dumusaue et al.ll2012l ; IGuedes et aTll2008l ). 

The AM amplitudes determined for the 19 systems at 
hand are well below 1 fias. This will put the systems in 
consi deration even beyond the reach of ESA 's GAIA mis- 
sion dHestroffer et al . 2010). However, recentlv lMalbet et alJ 
l|201ll ) proposed the Nearby Earth Astrometric Telescope 
(NEAT) which would be capable of resolving astrometric 
motion down to 0.05 /xas at a one a accuracy level. This in- 
strument would be able to identify habitable planets in most 
of the presented binary star systems. Such a mission would 
indeed be valuable, since AM does not only favor planet de- 
tection in binary configurations with Sun-like components - 
their HZs are further away form their host stars thus produc- 
ing larger AM amplitudes - it would more importantly grant 
observational access to all the planet's orbital parameters. 
Especially mutual inclinations are of interest in this C£LS6, ELS 
they could provide answers to many important problems re- 
garding planet-formation as well as migration in binary star 
syste ms (|Batygin et al.ll201ll ; [ThebaultJl201ll ; IWu fc Murray! 
l2003l) . 



1 In this case rms values have not only been time averaged, but 
they were also averaged over the planet's argument of pericenter. 
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6 POTENTIALLY TRANSITING SYSTEMS 

With an inclination of I ~ 90° with respect to the plane of 
the sky the systems HIP 14669, 31711, 64241 & 64797 could 
harbor transiting planets that still would be compatible with 
our assumptions of a planar binary planet configuration. As- 
suming non grazing transits, i.e. transits where less than the 
full planetary disc obscures the stellar surface during tran- 
sit, and neglecting entry as well as limb darkening effects, 
we can estimate the relative transit depth (TD) that the 
planet will cause in its host-star's photometric signal: 



TD ~ — 



(l) 



Hereby, R p and R+ denote the planetary and stellar radii 
respectively. Table [S] shows the relative transit depths for 
Earth-like planets in systems allowing for transits while still 
being close to planar. Even-though some stellar components 
are on the verge of being too bright to be observed by Kepler, 
the spacec raft's current perform ance (combined noise level 
« 29 ppm. lGilliland et"aD (|201ir 0. would allow for an Earth- 
like planet in circumstellar HZs to be found in all of these 
systems given sufficient observation time. 



which give inaccurate results for low inclinati on configura- 
tions such as di scussed in the s t udy at hand ( Lee fc Pealel 
2003). Similar to lEggletonl (|2006l ). ICorreia et a.1.1 (|201ll7 show 
that tidal interactions in inclined hierarchical triple systems 
can produce many different outcomes, especially when the 
component's changes in obliquity are taken into account. 
Their system's final states included transformations of ret- 
rograde to prograde motion and vice versa, a decay of mutual 
inclination and rapid circularization of the inner planetary 
orbit as well as tidally induced migration. As more detailed 
tidal interaction models require know ledge of the stellar radii 
|Correia et al.ll201lj ; lEggletorJl2006l l . the model dependence 
of r adii for M-dwarfs adds another source of uncertainty, see 
e.g. iMuirhead et~ai1 (|2012T i. 

The lack of accurate analytical tools to study the in- 
fluence of tidal interactions in planar S-Type configurations 
as well as the wealth of possible final states depending on 
the system's initial conditions put a detailed analysis of the 
planet-binary system's tidal evolution beyond the scope of 
this work. 



7 TIDAL LOCKING 

An orbital state, where the planet rotates around its own 
axis with the same speed as it orbits its host-star - much 
like the Moon around the Earth - is called 1:1 spin-orbit 
resonance. A star-planet system might evolve into such a 
state due to tidal interactions (see e.g. iMurrav fc DermottJ 
(1999)). Therefore, this state is often referred to as tidal lock. 
Since a tidal locking p otentially adds ad ditional instabilities 
to a planet's climate (|Kite et al.ll201l| ), regions where 1:1 
spin orbit resonance s occur are usually excluded fro m HZs. 
IKasting et all (| 19931 ) used an equation dating back to lPeald 
1 1977f ) to calculate the distance up to which a planet would 
be tidally locked in a time-span equal t o age of the Solar Sys - 
tem. Inserting such values as chosen in IKasting et all (|l993l ) 
the simple estimate reads: 

with ttl denoting the tidal locking radius in AU and m* 
the mass of the host star in Mq. Applying this estimate 
to our selected systems indicates that all HZs in M-M bi- 
naries fall at least partly in the Tidal Locking Zone. How- 
ever, tidal evolution of a planet in a binary system is much 
more involved than simple two body dynamics can account 
for, as the angular momentum transfer between the host- 
star-planet s ystem and the se condary needs to be included 
in the model. lEggletor] |2006) provides analytical estimates 
for the tidal evolution of stellar hierarchical triple systems 
showing that in fact many possible resonant states other 
than 1:1 spin-orbit locking exists for the inner pa i r alth ough 
with different degrees of stab ility. IWu fc Murravl (|2003t ) and 
iFabrvckv fc Tremainel (|2007T l investigated the possibility for 
tidal migration of planets d ue to mutual ly inclined massive 
pe rturbers via Kozai cy cles (|Kozailll962l ). Yet, as point out 
by ICorreia et al.l |201ll 'l , only quadrupolar secular expan- 
sions had been used to evaluate the planet's eccentricity, 



8 SUMMARY 

Apply ing the analytic methods presented in (Egg l et al.l 
l2012al lbT) we have shown that 17 out of 19 binary star sys- 
tems with well determined stellar and orbital parameters 
close to the Solar System allow for dynamically stable Earth- 
like planets in circumstellar habitable zones (HZs). Four of 
these habitable systems feature F, three feature G, six fea- 
ture K, and nine feature M class stars. Not surprisingly, 
M-M binary constellations offer the best chances for detect- 
ing planets in HZs via radial velocity observations. However, 
determining habitability in M star doublets may require ad- 
ditional con siderations such as t olerable stellar X-ray and 
EUV fluxes |Lammer et al.lfeoill ) or the system's potential 
for tidally locking the planet to its host star (see section [7]). 
Habitable planets in systems featuring G-type stars have 
RV amplit udes comparable to the ones found for a Cen- 
tauri AB (|Eggl et al.ll2012al : iGuedes et al.ll200Sl ). The sys- 
tems HIP 14699, 30920, 106972, 114922 or 116132 would be 
promising candidates to search for terrestrial planets in their 
HZs, as they offer best case RV semi-amplitu des compara- 
ble to a Centauri B b (|Dumusque et af]|2012h . Four of the 
17 systems would allow for transiting planets in HZs, which 
could be detected using current technology. Their mid tran- 
sit depths were estimated to lie between 44 and 369 ppm 
with planetary periods ranging from 235 to 476 days. As- 
trometric signal amplitudes for Earth-like planets in all the 
investigated systems' HZs are, in contrast, wel l below 1 /j,as. 
Ther efore, dedicated missions such as NEAT (Malbe t~et al.l 
l201lT ) will be required in order to detect habitable worlds 
in binary stars via astrometry. A sample of 19 systems does 
not offer the possibility to construct a reasonable statistical 
analysis on the number of potentially habitable binary star 
systems in the Solar neighborhood. More precise data on 
spectral types and orbital elements of nearby double stars 
are required in this respect. Nevertheless, our findings indi- 
cate that including binary star systems with 1 < ai < 100 
in observational campaigns has the potential to enhance our 
chances of finding habitable worlds. 
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Table 2. Critical semimajor axis (a cr n [AU], col. 3) for orbital stability and borders for the HZs ([AU], cols. 5-9) as 
measured for the respective host stars A&B are given for 19 binary star systems in the Solar neighborhood. Additionally, 
rms radial velocity (RV [cm/s]) and astrometric (AM [^tas]) signatures of terrestrial planets have been evaluated at the HZ 
borders. The conditions required for a planet to be withing the Averaged (AHZ), Extended (EHZ) and Permanent (PHZ) 
Habitable Zone are discussed in section [3] Dashed fields (-) represent cases where a given HZ border lies beyond the critical 
semimajor axis a cr ; t . Planets there would be on dynamically unstable orbits. 



HIP ID 


comp. 




inner AHZ 


inner EHZ 


inner PHZ 


outer PHZ 


outer EHZ 


outer AHZ 










0.306 


0.308 


0.310 


0.590 


0.596 


0.604 


HZ 








22.02 


21.95 


21.88 


16.04 


15.96 


15.86 


max RV 




A (M2) 


2.287 


15.39 


15.34 


15.29 


11.09 


11.03 


10.96 


rms RV 








0.107 


0.108 


0.109 


0.209 


0.211 


0.214 


max AM 


14669 






0.076 


0.076 


0.077 


0.146 


0.147 


0.149 


rms AM 






0.162 


0.162 


0.162 


0.316 


0.318 


0.320 


HZ 








36.00 


36.00 


36.00 


25.94 


25.86 


25.78 


max RV 




B (M4) 


1.806 


25.30 


25.30 


25.30 


18.12 


18.06 


18.00 


rms RV 








0.081 


0.081 


0.081 


0.159 


0.160 


0.161 


max AM 








0.057 


0.057 


0.057 


0.112 


0.112 


0.113 


rms AM 








0.086 


0.086 


0.088 


0.162 


0.166 


0.168 


HZ 








52.37 


52.37 


51.80 


38.90 


38.47 


38.26 


max RV 




A (M4V) 


0.865 


36.24 


36.24 


35.83 


26.43 


26.11 


25.95 


rms RV 








0.291 


0.291 


0.298 


0.557 


0.571 


0.578 


max AM 


60920 






0.237 


0.237 


0.242 


0.445 


0.456 


0.462 


rms AM 






0.027 


0.027 


0.027 


0.051 


0.051 


0.051 


HZ 








157.79 


157.79 


157.79 


115.07 


115.07 


115.07 


max RV 




B (M5.5V) 


0.470 


110.84 


110.84 


110.84 


80.36 


80.36 


80.36 


rms RV 








0.255 


0.255 


0.255 


0.487 


0.487 


0.487 


max AM 








0.210 


0.210 


0.210 


0.400 


0.400 


0.400 


rms AM 








0.886 


0.894 


0.902 


1.694 


1.724 


1.756 


HZ 








9.63 


9.59 


9.55 


7.09 


7.03 


6.97 


max RV 




A (G2V) 


8.351 


6.68 


6.65 


6.62 


4.83 


4.79 


4.74 


rms RV 








0.125 


0.126 


0.127 


0.243 


0.247 


0.252 


max AM 


31711 






0.087 


0.088 


0.088 


0.166 


0.169 


0.172 


rms AM 






0.316 


0.318 


0.320 


0.614 


0.618 


0.622 


HZ 








21.36 


21.29 


21.23 


15.42 


15.37 


15.32 


max RV 




B (K7Ve) 


5.848 


15.00 


14.95 


14.90 


10.76 


10.72 


10.69 


rms RV 








0.079 


0.080 


0.080 


0.155 


0.156 


0.157 


max AM 








0.056 


0.056 


0.056 


0.108 


0.109 


0.109 


rms AM 








1.581 


1.619 


1.697 


2.686 


2.686 


2.686 


HZ 








4.80 


4.75 


4.66 


_ 


_ 


_ 


max RV 




A (F3V) 


2.686 


3.19 
0.221 


3.15 
0.226 


3.08 
0.238 


- 


- 


- 


rms RV 
max AM 








0.176 


0.181 


0.189 








rms AM 






0.710 


0.718 


0.734 


1.340 


1.418 


1.456 


HZ 








8.76 


8.71 


8.62 


6.59 


6.44 


6.38 


max RV 




B (KOV) 


1.967 


6.03 


6.00 


5.93 


4.39 


4.27 


4.21 


rms RV 








0.154 


0.156 


0.160 


0.300 


0.320 


0.329 


max AM 








0.127 


0.129 


0.132 


0.241 


0.255 


0.262 


rms AM 








0.263 


0.263 


0.263 


0.515 


0.515 


0.517 


HZ 








8.79 


8.79 


8.79 


6.30 


6.30 


6.28 


max RV 




A (MOV) 


17.932 


6.20 


6.20 


6.20 


4.43 


4.43 


4.43 


rms RV 








0.265 


0.265 


0.265 


0.520 


0.520 


0.522 


max AM 


45343 






0.256 


0.256 


0.256 


0.501 


0.501 


0.503 


rms AM 






0.252 


0.252 


0.254 


0.494 


0.496 


0.496 


HZ 








9.08 


9.08 


9.04 


6.50 


6.48 


6.48 


max RV 




B (MOV) 


17.698 


6.41 


6.41 


6.38 


4.58 


4.57 


4.57 


rms RV 








0.259 


0.259 


0.261 


0.509 


0.511 


0.511 


max AM 








0.250 


0.250 


0.252 


0.491 


0.493 


0.493 


rms AM 


51986 


A (F4IV) 


0.545 


HZ 


B (F3) 


0.448 


HZ 


58001 


A (AOVe) 


2.828 


HZ 






0.639 
10.31 


0.653 
10.21 


0.775 
9.46 


1.263 
7.85 


1.331 


1.331 


HZ 

max RV 




B (K2V) 


1.331 


6.98 
0.100 
0.080 


6.91 
0.102 
0.082 


6.34 
0.123 
0.097 


4.97 
0.210 
0.159 






rms RV 
max AM 
rms AM 
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Table 3. Continuation of Table [2] Radial velocity (RV) amplitudes are given in [cm/s], astrometric (AM) amplitudes in 
[/(as]. The critical planetary semimajor axis a cr n as well as the HZ borders are given in [AU]. 



HIP ID 


comp. 




inner AHZ 


inner EHZ 


inner PHZ 


outer PHZ 


outer EHZ 


outer AHZ 










1.354 


1.465 


1.465 


1.465 


1.465 


1.465 


HZ 








8.19 


_ 


_ 


_ 


_ 


_ 


max RV 




A (F5V) 


1.465 


4.82 
0.209 


_ 
_ 


_ 

_ 


_ 
_ 


_ 

_ 


_ 
_ 


rms RV 
max AM 


64241 






0.126 


_ 


_ 






_ 


rms AM 






1.002 
9.76 


1.056 
9.58 


1.226 
9.15 


1.339 
- 


1.339 
- 


1.339 
- 


HZ 

max RV 




B (F6V) 


1.339 


6.05 
0.173 
0.109 


5.90 
0.184 
0.115 


5.47 
0.219 
0.133 


- 


_ 

_ 


_ 
_ 


rms RV 
max AM 
rms AM 








0.472 


0.472 


0.474 


. 


0.926 


0.926 


HZ 








15.48 


15.48 


15.45 


11.08 


11.07 


11.07 


max RV 




A (K1V) 


23.212 


10.93 


10.93 


10.91 


7.81 


7.80 


7.80 


rms RV 








0.179 


0.179 


0.180 


0.351 


0.352 


0.352 


max AM 


64797 






0.126 


0.126 


0.127 


0.248 


0.248 


0.248 


rms AM 






0.263 


0.263 


0.263 


0.517 


0.517 


0.517 


HZ 








24.51 


24.51 


24.51 


17.50 


17.50 


17.50 


max RV 




B (M1V) 


18.564 


17.32 


17.32 


17.32 


12.35 


12.35 


12.35 


rms RV 








0.140 


0.140 


0.140 


0.275 


0.275 


0.275 


max AM 








0.099 


0.099 


0.099 


0.194 


0.194 


0.194 


rms AM 








0.339 


0.341 


0.345 


0.645 


0.655 


0.667 


HZ 








12.20 


12.17 


12.10 


8.99 


8.92 


8.85 


max RV 




A (M0.5) 


4.289 


8.48 


8.46 


8.41 


6.15 


6.10 


6.05 


rms RV 








0.081 


0.081 


0.082 


0.156 


0.158 


0.161 


max AM 


66492 






0.072 


0.072 


0.073 


0.137 


0.139 


0.142 


rms AM 






0.227 


0.229 


0.231 


0.439 


0.443 


0.449 


HZ 








16.40 


16.33 


16.26 


11.92 


11.87 


11.80 


max RV 




B (Ml. 5) 


3.835 


11.46 


11.41 


11.36 


8.24 


8.21 


8.15 


rms RV 








0.066 


0.066 


0.067 


0.129 


0.130 


0.132 


max AM 








0.059 


0.060 


0.060 


0.114 


0.115 


0.117 


rms AM 








0.486 


0.490 


0.496 


0.916 


0.934 


0.952 


HZ 








11.47 


11.43 


11.36 


8.51 


8.43 


8.36 


max RV 




A (K4V) 


4.503 


7.95 


7.92 


7.87 


5.79 


5.74 


5.68 


rms RV 








0.155 


0.157 


0.159 


0.298 


0.304 


0.310 


max AM 


67422 






0.130 


0.131 


0.133 


0.245 


0.250 


0.255 


rms AM 






0.398 


0.400 


0.404 


0.754 


0.766 


0.780 


HZ 








13.37 


13.34 


13.27 


9.85 


9.78 


9.70 


max RV 




B (K6V) 


4.212 


9.30 


9.27 


9.23 


6.75 


6.70 


6.64 


rms RV 








0.142 


0.143 


0.144 


0.273 


0.277 


0.282 


max AM 








0.119 


0.120 


0.121 


0.226 


0.229 


0.234 


rms AM 


84425 


A (F7V) 


1.024 


HZ 






0.635 
12.54 


0.667 
12.32 


0.797 
11.67 


0.835 


0.835 


0.835 


HZ 

max RV 




B (G2V) 


0.835 


7.82 
0.082 
0.056 


7.63 
0.087 
0.059 


6.98 
0.107 
0.071 








rms RV 
max AM 
rms AM 








0.535 


0.543 


0.551 


1.003 


1.029 


1.057 


HZ 








8.34 


8.29 


8.23 


6.25 


6.17 


6.10 


max RV 




A (G8V) 


4.275 


5.73 


5.69 


5.65 


4.19 


4.13 


4.08 


rms RV 








0.240 


0.244 


0.248 


0.461 


0.474 


0.487 


max AM 


84720 






0.213 


0.216 


0.219 


0.399 


0.410 


0.421 


rms AM 






0.242 


0.242 


0.244 


0.462 


0.468 


0.474 


HZ 








15.40 


15.40 


15.33 


11.27 


11.20 


11.14 


max RV 




B (MOV) 


3.364 


10.75 


10.75 


10.70 


7.78 


7.73 


7.68 


rms RV 








0.170 


0.170 


0.172 


0.329 


0.333 


0.337 


max AM 








0.153 


0.153 


0.154 


0.292 


0.296 


0.300 


rms AM 


87895 


A (G2V) 


0.371 


HZ 


B (K2V) 


0.312 


HZ 
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Table 4. Continuation of Table [2] Radial velocity (RV) amplitudes are given in [cm/s], astrometric (AM) amplitudes in 
[//as]. The critical planetary semimajor axis a cr u as well as the HZ borders are given in [AU]. 



HIP ID 


comp. 


&cvit 


inner AHZ 


inner EHZ 


inner PHZ 


outer PHZ 


outer EHZ 


outer AHZ 










1.289 


1.311 


1.337 


2.421 


2.505 


2.581 


HZ 








3.71 


3.68 


3.65 


2.79 


2.74 


2.71 


max RV 




A (F8V) 


5.623 


2.54 


2.51 


2.49 


1.85 


1.82 


1.79 


rms RV 








0.185 


0.189 


0.193 


0.358 


0.371 


0.383 


max AM 


93825 






0.168 


0.170 


0.174 


0.315 


0.326 


0.336 


rms AM 






1.254 


1.274 


1.300 


2.358 


2.438 


2.512 


HZ 








3.79 


3.76 


3.72 


2.84 


2.80 


2.76 


max RV 




B (F8V) 


5.575 


2.59 


2.57 


2.54 


1.89 


1.86 


1.83 


rms RV 








0.183 


0.186 


0.190 


0.353 


0.365 


0.377 


max AM 








0.165 


0.168 


0.171 


0.311 


0.321 


0.331 


rms AM 


101916 


A (GUV) 


0.754 


HZ 






0.155 
38.16 


0.159 
37.73 


0.233 
31.89 


0.261 
30.40 


0.381 


0.381 


HZ 

max RV 




B (K2IV) 


0.381 


25.63 
0.045 
0.031 


25.31 
0.046 
0.032 


20.91 
0.069 
0.047 


19.75 
0.078 
0.053 


_ 


_ 


rms RV 
max AM 
rms AM 








0.322 


0.330 


0.338 


0.588 


0.616 


0.640 


HZ 








20.73 


20.50 


20.28 


15.94 


15.65 


15.41 


max RV 




A (M2) 


1.034 


13.99 


13.82 


13.65 


10.35 


10.12 


9.93 


rms RV 








0.074 


0.076 


0.077 


0.139 


0.147 


0.153 


max AM 


106972 






0.053 


0.054 


0.055 


0.096 


0.101 


0.105 


rms AM 






0.179 


0.183 


0.185 


0.339 


0.351 


0.359 


HZ 








31.41 


31.08 


30.92 


23.35 


22.99 


22.76 


max RV 




B (M4) 


0.865 


21.63 


21.40 


21.28 


15.73 


15.46 


15.28 


rms RV 








0.053 


0.055 


0.055 


0.103 


0.107 


0.110 


max AM 








0.039 


0.040 


0.040 


0.074 


0.077 


0.078 


rms AM 








0.239 


0.245 


0.251 


0.435 


0.455 


0.473 


HZ 








24.59 


24.31 


24.05 


18.89 


18.53 


18.24 


max RV 




A (Ml) 


0.897 


16.60 


16.40 


16.20 


12.30 


12.03 


11.80 


rms RV 








0.050 


0.051 


0.053 


0.094 


0.099 


0.103 


max AM 


114922 






0.037 


0.038 


0.039 


0.068 


0.071 


0.074 


rms AM 






0.266 


0.272 


0.282 


0.480 


0.504 


0.528 


HZ 








22.83 


22.60 


22.23 


17.67 


17.32 


17.00 


max RV 




B (M2) 


0.924 


15.33 


15.16 


14.89 


11.41 


11.14 


10.88 


rms RV 








0.053 


0.055 


0.057 


0.100 


0.105 


0.111 


max AM 








0.039 


0.040 


0.042 


0.071 


0.075 


0.078 


rms AM 








0.158 


0.158 


0.158 


0.310 


0.310 


0.312 


HZ 








30.77 


30.77 


30.77 


22.02 


22.02 


21.95 


max RV 




A (M4) 


10.619 


21.72 


21.72 


21.72 


15.51 


15.51 


15.46 


rms RV 








0.204 


0.204 


0.204 


0.401 


0.401 


0.404 


max AM 


116132 






0.165 


0.165 


0.165 


0.323 


0.323 


0.325 


rms AM 






0.076 


0.076 


0.076 


0.152 


0.152 


0.152 


HZ 








60.97 


60.97 


60.97 


43.14 


43.14 


43.14 


max RV 




B (M5) 


7.091 


43.07 


43.07 


43.07 


30.45 


30.45 


30.45 


rms RV 








0.184 


0.184 


0.184 


0.369 


0.369 


0.369 


max AM 








0.149 


0.149 


0.149 


0.298 


0.298 


0.298 


rms AM 



Table 5. Transit depths (TO), visual brightness (V, WDS) 
and planetary period (P p ) ranges are given for potentially 
transiting planets in the HZs of those selected binary systems 
with I « 90° . 



HIP ID 


comp. 


V [mag] 


TD [ppm] 


P P [D] 


14669 


A 


10.32 


128 


270.80 - 380.35 


B 


12.5 


369 


235.67 - 331.14 


31711 


A 


6.32 


78 


338.38 - 476.35 


B 


9.84 


187 


270.91 - 380.10 


64241 


A 


4.85 


44 


372.23 


B 


5.53 


79 


346.09 - 382.82 


64797 


A 


6.66 


171 


294.45 - 412.52 


B 


9.5 


198 


260.18 - 364.81 
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Figure 1. Habitable zones of 19 S-Type binary star systems in the Solar neighborhood are shown. The light gray regions (blue online) 
denote zones of permanent habitability (PHZ), medium gray (green online) Extended (EHZ) and dark gray (yellow online) Averaged 
Habitable Zones (AHZ), see section [3] Black (red online) are regions where the planet either rec eives too much, or too litt le radiation 
to keep atmospheric temperatures stable. The striped areas are zones of dynamical instability jHolman fc Wiegerd fl999h . Left: HZs 
around the system's p rimary star are shown (S-Type A), right: habitability of regions around the secondary are investigated (S-Type B) 
(Whitmirc et al . 1998). The dashed T symbols i ndicate the inner, the full symbols the outer border of the classical HZ as defined by 
iKasting et al. I (11993!) &: lUnderwood et alj d2003l '). In most cases, the AHZ and the classical HZ coincide well as was pointed out in 
(Eggl et al. 2012b), except for the systems HIP 58001 and 101916 where the considerable luminosity of the brighter companions shift 
the HZs of the S-Type B configurations to larger planetary semimajor axes. Evidently, 17 out of the 19 investigated systems allow for 
dynamically stable terrestrial planets within HZs around at least one of its binary's components. 



APPENDIX A: MAXIMUM AND RMS SIGNAL AMPLITUDES 

Following iBeauge et al.l l|2007l ) & lEggl et all (|2012al ) the RV amplitude a planet causes on its host star is given by 

\/Gmi sin I e coscj + cos(/ + uj) 



RV 



\Jmv + rm 



(Al) 



where G denotes the gravitational constant, m,n and mi are the host-star's and planet's masses. The quantities a, e, I and uu 
denote the planet's semimajor axis and eccenticity, the system's inclination to the plane of the sky and the planet's argument 
of pericenter respectively. The planet's true anomaly is represented by /. We can write the maximum possible radial velocity 
(RV) amplitude caused by a terrestrial planet in a circumstellar orbit around one binary component as follows: 



max RV = 



t/Q mi sin I (1 + e T ' 



\/mo + mi ^Ja(l - ( e max ) 2 ) 



(A2) 



The maximum possible eccentricity the planet can acquire due to gravitational interaction with the double star is denoted by 



l2012al ) 



(Egg l et al.ll2012bl ). Expressions for the root mean square (rms) values of the RV signal are given as follows l|Eggl et alj 



rms RV 



1 

2^ 



I 1/2 



RV'dMdu 



\[Q mi| sin/| 
*yj2a(mQ + mi 



(A3) 
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In a similar manner we can use the formalism applied in iPourbabd |2002l ) to determine maximum astrometric (AM) signal 
strengths: 

„„, fia(l + e max ) , tJ . 

max AM = — - — - (A4) 

where d is the observers distance to the observed system, and /j, — mi/ (mo + mi) is the planet-host-star system's mass ratio. 
The astrometric rms amplitudes are given by: 



rms AM = ^~ 
2d 



3+?<e 2 ) + (l + ^e^))cos(2/) 



1 1/2 

(A5) 



Here, (e 2 ) is the averaged squared planetary eccentricity. An analytic expression for (e 2 ) can be found in iGeorgakarakol (|2003l . 
2005). For a detailed derivation of all rms and maximum signal amplitudes see lEggl et al. | (|2012al ). 



